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This paper describes a time-sharir.g FORTRAN program whi^h allowed a student to design 
and perform a variety of classical thermodynamic experiments on three common models of a 
working substance. An important goal of this project was to develop a simulation which gave 
th^- student maximum freedom and responsibility in the design of the experiment. and provided 
only v he primary experimental results. In this way the student obtained valuable practice in 
planning and data analysis. The particular purpose of the project was to develop a program 
which the student could use to perform simulated thermodynamics experiments. The decisions 
about experimental variables to -Measure and control, and the set of v. lues to use, must be 
made by the student; the data analysis is left to the student, tven if it could easily have 
been aided to the program. An attempt has been made to provide, as output, only those 
thermodynamic variables which are directly measurable in the lahoiat-ory. As an example, an 
assignment may be to determine the specific heat vs. temperature curv P for the substance 
under study. Specific heat is not directly measurable and it is the student's responsibility 
to determine whac parameters to measure, what conditions to us->, and what analysis to mak^. 

EEM^s m_D escr ipt i on * 

PVTEXP is a program, written in Call-A-Computer FORTPAN, which simulates three of *h^- 
most popular models for a real substance. These are (a) the ideal gas, (b) Van der Waal's 
gas, and (c) a solid with a Debye specific hea 1 - curve. The program is used at our 
interactive, timesharing terminal, but this mode of operation is not felt to be necessary for 
successful operation. 

The "User In str uct ions M attached to this paper are those given to our students. The 
instructions purposely omit the names of the models and refer to them only by number. This 
is done so that familiarity with "expected results" will not prevent discovery. Following 
thase instructions, the student selects a working substance by number; chooses a fixed 
parameter (i.e., an isothermal, isobaric, isochoric or adiabatic process) with specified 
value; chooses to vary one of the remaining parameters; and enters an initial value, final 
value, and measurement interval. If he has chosen heat flow as one parameter, he will be 
asked to enter an additional P-V-T coordinate of the initial state since two such coordinates 
are req?:i-ed to specify a point on the state surface. For each of the models, no^ counting 
differences due to numerical values, there are fifteen ways of specifying the desired 
experiment. After appropriate selection of values and subsequent analysis of results, this 
set of 45 experimental choices may be combined in a wide variety of ways and made to answer 
many types of questions. The printed results are pressure in atmospheres, volume in liters, 
temperature in °K, and cumulative heat flow in calories at each step of the requested 
process. The units may be easily changed by an appropriate change of conversion factors in 
the DATA lines. In addition, the user will find it an easy matter to change model parameters 
(such as the Van der waals constants or the Debye temperature) and study the effects of these 
changes. 

furthermore, since the Program is essentially a series of subroutines, it is possible to 
introduce other models with relative ease. It does require, however, an eguation of state 
(with algorithms for finding, P, V and T) and an ex-pres-sion. for Cp or C v in terms of these 
state variables. However, if no convenient algorithm for finding rV-(P,T), for example, is 
available, then a new model may still be introduced with the restriction that V is +o be 
specified by the experimenter, either as a fixed or stepped parameter. 



m_A ££l icat i on 

At Wofford, the program has found application at two places in the curriculum. One of 
these is with a class of freshman "non-scientists" who were given very little background 
infornation on thermodynamics. They had done some (non-thermo) laboratory experiments, 
statistical analysis of data, curve fitting, and a small amount of computer programming. The 
main theme in this course was the quantitative ' treatment of data in science. The 
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t hemodynamics program describee 1 here was introduced as an example of computer simulations of 
nature. Applications in this course included: 

1. After a brief lecture on the nature of an equation of state and how such equations 
may be obtained experimentally, these students were asked to design experiments on 

?' "substance number 1 " to obtain some information about its* equation of state. Some 

"discovered" Boyle's law, some found Charles 1 law, and a few found combined results 

and reported something like T - CPV with a value of C and some indication of its 
variabili t y. 

2. After a brief discussion of specific heat and the fact that there are different 
specific heats for different processes, the students wer=> asked to find, for 
substance number 1, C v at an assigned temperature. 

3. Aft^r a classroom discussion of thermal expansion, the qroup was asked to select 
and perform experiments to determine the pressure dependence of the coefficient of 
expansion tor substance number 2. Each student took data at one or more values of 
pressure ar.tl a riot was made of the group results. Subsequent discussion of the 
variability in these results proved fruitful for those who stuck with this problem. 

I olrin to make similar uses cf the program with a new group of students in the spring 
semester . 

The second application is in the thermodynamics course. A typical physics course in 
+ hermod ynamics has no lab, but the experimental procedures for determination of specific 
heat, enthalpy, and other functions are important research methods. furthermore, students 
are asked to work with several new abstract functions without the experience of seeing in a 
personal way how such functions are obtained frorr experiments. Therefor* 3 , this course is an 
ideal one for the introduction of simulated experimental planning and analysis. However, a 
literature search and other inquiries revealed very little work in the area and no programs 
with ^he desired capability. 

flith PVT^TXP, students can be asked to observe system behavior, evaluate parameters, 
determine functional relationships, or investigate properties of th s rmod y namic cycles. A 
list of example assignments is given below. At Wofford, our experience with advanced 
students using, the program is limited by the number cf students who take the course. We have 
not exhausted the possibilities of the present program, but already recognize the need for an 
expanded version which will simulate phase changes and still allow the freedom of 
experimental design available with the present version. 

5§.IH£l§_P£9.^l£!!l_Ii§ss_Using_P V TEXP 

1. Choose substance number 1 and evaluate beta and kappa near STP. 

2. Choose substance number 2 and plot kappa as a function of temperature from t 50°K to 
1500°K at P = 1 a tin. 

3. Choose substance number 3, measure and plot bexa as a function of temperature from 
50°K to 2000°K. 

U. Experimentally determine values for the Van der Waals constants, a and b, using 
substance number 2. 

5 . . Choose substance number 2 and plot PV/PT ageinst P from 1 to 800 atm. at several 

temperatures. 

6. Find the work done in isothermally changing one mcle of substance number 1 at 300°K 
and 1 atm. to 300°K and 5 atm. " ' 

7. Determine C p for substance number 1 in the neighborhood of STP. 

3. Determine and plot C p of substance number 3 as a function of temperature from 10°K 
to 400°K. 

9. Determine the enthalpy difference between 1 mole of substance number 2 at 300°K and 

at 50Q°K, both at 2 atm. 
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13. Determine and plot tf as a function of T for substance number 1 from T = 100°K to T * 
= 200°K at 1 atm. (Use an arbitrary reference level.) 



11. Use substance number 1 and substance number 2 in separate throttling processes. 
Determine the Joule-Thompson coefficient of each substance in a region of interest 
to you. 

12. Use substance number 1 in a Carnot cycle cf your choice and calculate th« 
efficiency from W and Q. Compare your result with the theoretical value. 

13. Choose two Carnot cycles which have one boundary in common. With substance number 
2 as the working fluid,, experimentally determine the efficiency of each cyle and of 
the two cycles together. 

1^. -Choose a Carnot cycle and an Otto cycle with the same temperature limits. 
Determine and compare their efficiencies with substance number 1 as the operating 
fluid. 



15. Experimentally determine the Debye temperature for substance number 3. 

16. Plot entropy as a function of temperature for substance number 3 from 0°K to 200°K. 



PVTEXP 



User Instructions 



PVTEXP is a FORTRAN program to allow the user to design experiments and run selected 
processes on the P-V-T surfaces of several substances. The units for all input and output 
quantities are: pressure in atmospheres, volume in liters, temperature in degrees Kelvin, and 
heat flow in gram-calories. One gram-mole of substance is represented in each ci.se. A 
process on the P-V-T surface may be chosen to be: 

1. constant pressure (isobaric) 

2. constant volume (isoboric) 

3. constant temperature (isothermal) 
constant Q (a diaba tic) 

NOTE: q is not a state variable in thermodynamics. Only changes in Q have meaning. 
However, in this Drogram, an adiabatic change is specified as Q = 0 and a process of 
controlled heat input may be specified as (for example) Q = 0, 200, 10 which is taken to 
mean "add heat in steps of 'l 0 calories until 200 calories have been added." 

There are three substances in the present program. The user selects one by responding 
1, 2, or 3 to the first question which is "WHICH SUBSTANCE — TYPE NUMBER? " . 

In response to the question "SPECIFY FIXED PARAMETER?", a constant value for one of the 
four parameters is entered by typing an equation of the form T = 273. You may use any one of 
the variables P, V, T, or Q and assign its constant value in this way. (Other examples: P = 
1.2, V = 2U.3, Q = 0.) In ent3ring equations, do not use a space anywhere before the equal 
sign . 

The next question is "SPECIFY STEPPED PARAMETER WITH INITIAL AND FINAL VALUES AND STEP- 
SIZE?". A proper reply has the form P = 1,2,. 2 which specifies a. process beginning at P = 1 
atmosphere,, moving along the P-V-T surface (in the direction specified earlier) until P = 2 
atmospheres, printing values of P,V, t and heat flow at each step of 0.2 atmospheres. Any 
one of the four variables may be specified in this way. ^ 

Two coordinates are required to locate the point from which you wish to start your 
process and Q is not a coordinate. If Q is specified as either the fixed or the stepped 
parameter, the next question will be "SPECIFY ANOTHER P-V-T COORDINATE OF T HE STARTING 
POINT?" . This is answered with another equation like V 25. 

Four replies to these questions have specified a path on the P-V-T surface — the starting 
point, the "direction", and the svopping criterion. Pressure, volume, temperature and net 
heat flow are then listed for points along this pa^h, using the interval specified as "stpp- 
size" earlier. 

When the tabulation is complete, the computer returns to the question "WHICH SUBSTANCE — 
TYPE NUMBER?" and you may continue with a., new problem or you pay wish to enter the previous 
final values of P r V, and T to begin a new part of a process under study. In +his way, you 
may study closed cycles. Enter a zero (or an S) to stop the program. 
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Sample Run 
PVTEXP 02H7 SAT. 03-24-73 



WHICH SUBSTANCE — TYPE NUMBER?2 
SPECIFY FIXED PARAMETER?!* 200 

SPECIFY STEPPED PARAMETER WITH INITIAL AND FINAL VALUES AND STEP- SIZE 



?P*. 5* 1. 5# .25 

PRESSURE VOLUME TEMPERATURE HEAT FL0W 

< ATM. ) (LITERS) < DEG. KELVIN) <CAL0RIES) 

.50 32.8356 200.00 0*00 

.75 21.3943 200.00 -161.3181 

1.00 16.4237 200.00 -275.7905 

01.25 1 3.1413 200.00 - 3 64.58 73 

1.50 10.9 53 200.00 -437.1387 



WHICH SUBSTANCE — TYPE NUMBER? 1 
SPECIFY FIXED PARAMETER? V* 25 

SPECIFY STEPPED PARAMETER W3 TH INITIAL AND FINAL VALUES AND STEP- SIZE 
?T* 200* 300,20 



PRESSURE 


VOLUME 


TEMPERATURE HEAT FLOW 


< ATM. ) 


(LITERS) 


<DE6. 


KELVIN) <CAL0RIES) 


.6565 


25.00 


200.00 


o.oo 


.7221 


25.00 


220.00 


59 . 6 478 


.7878 


25-00 


240.00 


119.29 55 


.3 534 


25.00 


260.00 


178.9 433 


• 9 I? 1 


25.00 


280.00 


238.5911 \ 


.98 47 


25.00 


300.00 


298 . 2388 



WHICH SUBSTANCE — TYPE NfFiBER?3 
SPECIFY FIXED PARAMETER7P* I 

SPECIFY STEPPED PARAMETER WX TH INITIAL AND FINAL VALUES AND STEP-SIZE 
?T=» 100*1000,200 



PRESSURE 


V0LUME 


TEMPERATURE HEAT FLOW 


< ATM.) 


(LITERS) 


<DE6. 


KELVIN) (CALORIES) 


1.00 


• 007 


100. 00 


0.00 


1.00 


.0071 


300.00 


1127.7404 


1.00 


• 0072 


500.00 


2347.7648 


1.00 


.0073 


700.00 


3609.0076 


1.00 


.0073 


900.00 


4910.6003 


1.00 


.007 4 


1100.00 


6256.9703 



WHICH SUBSTANCE — TYPE NUMBER70 



AT LINE N0. 210: STOP » RAN 216/6 SEC. 



Program Listing 



PVTEXP 



100C VERSION OF 22 AUGUST 1972* DAN W. OLDS* W0FF0RD COLLEGE 
110 DIMENSION S<4)»WC9) 
120 REAL KAPPA 

130 EQUIVALENCE < S< 1 )* P) * < S< 2) * V)* ( S< 3) * T)* C S< 4) * Q) 

1 40 \COMM0N FF< 4) 

150 1F0RMATC 1A1*1X*I15)J 2F0RMATC 1 Al* IX) 

160 READ* W» FF 

170 R-8314.9 

180 AS»248 00*B5=.0266 

190 BEGINt PRINT* t » 

200 PRINT W WHICH SUBSTANCE — TYPE NUMBER*** I INPUT* N 1 
210 I F< -Ml) OK J STOP 
220 LETTR(M): 

230 LETTR»* 32-M0DCM/4O9 6*20) )/7l RETURN 
240 V0LCP*T)» 

2 50 60T0C 50* 60* 70)N1 
260 5OV0L a R*T/PJ RETURN 
270 60 V5*R*T/P* ITER»0 

280 61 V6=R*T/P+B5-A5/<P*V5)*A5*B5/<P*V5*V5) 
290 IFClE-5-ABSf 1 V6/V5) > 62J VOL* V6J RETURN 
300 621 F( 2 5- ITER) 6 31 1 TER= I TER+ 1 J V5=V6*G0T061 

310 70V0L*.Q071*EXP<-l.348 43E-2*4. 1 5003E- 5*T+ 1 . 15E-8*T*T-P* 
320 +< 6. 7 111 E- 1 2+2. 7 1 405E- 1 5* T) ) 3 RETURN 

330 63PRINT* "VOLUME SOLUTION DOES NOT CON VERGE-- A VOI D THIS REGION OR 

3 40 ♦ SPECIFY VOLUME. "I STOP 
350 PRES<V,T)I 

360 GOTOC 150* 160* 1 70)N1 

370 150 PRES*R*T/ VI RETURN 

380 160 PRES=R*T/< V-85)-A5/< V*V)J RETURN 

390 170PRES=< AL0GC V/.007D+ 1. 34843E-2-4. 1 5003E- 5*T- 1 . 1 5E-8* T*T) / 

400 +C-6.711 1E-12+2.71405E-15*T> 

410 TEM(P,V)« 

420 GOTOC 250* 260* 270) N 1 

430 250 TEM=»P*V/RJ RETURN 

440 260 TEM*<P+A5/Vt2)*< V-B 5) /Rl RETURN 

450 270A1»4. 1 5003E- 5-2. 71 405E- 1 5*Pl AO*- 1. 348 43 E- 2- 6. 71 1 1E-12*P 

4cA0 TEM=<-AH-S9RTF<Al*Al-4.6E-8*<A0-AL0G< V/.0071) ) ) ) /2. 3E-8 

470 CBK(P*V*T)I 

480 G0T0 ( 440* 460* 480)N 1 

490 440 CV=1.5*R 

500 BETA* 1 /T 

510 KAPPAM/P 

520 CP*C^R 

530 RETURN \ 
540 4601 F< 50- T) 462J CP=3. 5* RJ GO TO 464 

5 50 462U» 3070. /TJ CP» R* < 3- 5* < 2* U/< EXP< U) - EXP( - U) ) ) » 2) 

560 464CP a CP+2*A5*P/< R* T* T) - < 6*A5*85/< R* 2* T» 3) - 6*A5t 2/C R» 3* T» 4) )*P» 2 

570 BETA=R/<P*V-A5/< V* 2* A5*B5/Vt 2) ) 

580 KAPPA = l/<-2*A5/< V*V) + R*T* V/C V-P"; >2) 

590 CV*CP-T*V*BETAt2/KAPPA 



PVTEXP CONTINUED 



600 RETURN 

610 480 BETA«41.5E-6+T*23E-9-P*.271405E-14 I 
620 KAPPA- 6. 7 1 1 1 E- ! 2+ 2. 7 1 405E- 1 5*T 

630 X«215./TlIF<X-l)482J DEB»O.J NT* 2* 10/XJ 8)0481/ CN«1#NT 

640 481DEB»DEBMX*3+3*Xt2/CN+6*X/CN»2*6/CNt3)*EXPt-X*CN)/CN 

650 DEB*3*<3« 1 41 59 26t 4/1 5-DEB) /Xt 01 G0T0483 

6 6 0 48BDEB»-3*X/8+l»+X*X/20.-Xt 4/1 680-Xt 6/90720 
67 0 483CV=»3*R*C4>i>DE8-3*X/CEXPCX)-l)) 

680 CP*CV*T*V*BETA»2/KAPPA 
690 SL0PECN2)! 

7 00 CALL CBK t P# V# T) 

7 10 G0T0CA#B#C#D*E*F)N2 

720 At SL0PE*CP/< BETA* V)f RETURN 

7 30 Bf SLOPE* CP* RETURN 

7 40 CiSL0PE"CWKAPPA/BE7AJ RETURN \ 

7 50 Of SLOPE* CVJ RETURN f 
760 El SL0PE»-T*V*BETAJ RETURN 

770 Pi SLOPE* 1*BETA/KAPPAl RETURN 

780 PVTCDl 

790 G0T0(P#V#T)L 

800 Pi PVT*PRESCV# T) J RETURN 

810 VI PVT*V0LCP#T) J RETURN 

820 Tl PV7>TEM(P# V)f RETURN 

830 FGCX#Y#Z)t 

8 40 N3«N2-6 

850 60T0<7>8#9# 10, 1 1, 12)N3 

8 60 7 CALL CeK<Z,Y,X) 
870 F*Y*X*BETA/CP 
880 G*-Y*CV*KAPPA/CP 
890 RETURN 

900 8 CALLCBKCY,Z,X) 
910 F«-BETA*X/(CV*KAPPA) 
920 G*-CP/<KAPPA*Z*CV) 

9 30 RETURN 

9 40 9 CALL CBK € Y, X# Z ) 

9 50 P«-KAPPA*CV/CBETA*Z)JG=CP/CZ*X*BETA) J RETURN 

960 10CALLCBKCP,Y,X) 

970 F* 1 . /CPI G*Y*BETA/CPI RETURN 

980 1 1CALLCBKCY,V,X)JF*1./CV*G=8ETA/CKAPPA*CV) J RETURN 

990 12CALLCBKC Y» X, T)l F*K APPA/C BETA*T) J G*- 1 . /< T*X*BETA) J RETURN 

1000 <JUTPCX,Y,Z)I 

1010 N3*NB-6 ?o 
1020 G0TCC2O, 30, 40* 45, 47, 49 ) N 3 
1030 20CALLL IN E<Z#Y,X, 0)1 RETURN 
1040 30CALLLINECY»Z,X,0) J RETURN 
1050 40CALLLINECY, X, Z#0)J RETURN 
1060 45CALLLINECP,Y,X,Z) J RETURN 
1070 47CALLLINEC Y, V, X,Z) J RETURN 
1080 49CALLLINEC Y, X, T, Z)J RETURN 
1090 END0 UTP 
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PVTEXP C0MTINUED 



1100 0Kt 

1110 PRINT" SPECIFY FIXED PARAMETER"* J IMPUT1* I* VALl I«LETTRC I > J SC I > =» VAL/FFC I > 
1120 PRINT M SPECIFY STEPPED PARAMETER WITH INITIAL AMD FINAL VALUES AND STB 
1 130 INPUT2*J*INPUT# VAL* ST9P* STEPJ J*L ETTRC J> J SC J> * VAL/FFC J> J ST0P»ST9P/FFC J> 
1140 STEP«STEP/FFCJ> 

1 150 I FCC 3. 5-I>*(3. 5-J>>70JN2=»3*I+J-4JN2*N2-N2/4J SC 6-I-J>»PVTC 6-1 -J> 

1 160 CALL TITLE 

1170 Q«0l CALLLINECP, V* T> 0> 

1180 VAR*START*SCJ> 

1 190 L»6-I-J 

1200 597DQ«0 

1210 DO 598* JJ»0*8 

1220 SCJ>» VAR+JJ*STEP/8 

1230 SCL>«PVTCL> 

1240 598 DG»SL0PECN2>*WCJJ+1> + DQ 
1250 0*Q+STEP*DQ/24 
1260 599 CALLL IN EC P* V* T* Q) 
1270 VAR«VAR+STEP 

1 280 IFC ABSC STOP* START) - ABSC VAR*. 001* STEP- START) > BEGIN* BEGIN* 59 7 
1290 70PRINT**SPECI FT ANOTHER P-V-T COORDINATE OF THE STARTING POINT** 
1300 INPUT£*K* VAL»K»LETTRCK) 

1 310 SCK>*VAL/FFCK>J CALLTITLEJ SC 10-J-K-I >»PVTC 10-J-K-I) 

1320 IFC3-J>80IN2=*J+6IX0*SCMAX0CK*6-K-J>>I Y0=SCMIN0CK* 6-K-J>>! 60T07 5 

1330 80N2«I+9IXO*SCMAXOCK*6-I -K) > J Y0» SCMINOCK* 6-I-X>> 

1 340 7 5ZO=SCJ>»CALL0UTPCXO*YO,ZO>lH=STEP 

1350 CALLFGCXO*YO»ZO> 

1360 XX»XO+H*F 

1370 YY«YO+H*G 

1330 Z»ZO*H 

1390 CALLFGCXX*YY*Z> 

1400 Xl=. 5*CX0+XX*H*F> 

1410 Y1«.5*CY0+YY+H*G> 

1420 CALL0UTPCX1*Y1*Z> 

1430 55CALLFGCX1,Y1,Z> 

1440 X20*X2*XO+2.*H*F 

1 450 Y20»Y2*Y0+2*H*G 

1460 F1»F* Gl a G 

1470 10CALLFGCX2,Y2*Z+H> 

1480 X2NEW*X1*.5*H*CF1 + F> 

1490 Y2NEW»Yl + » 5*H*CG1 + G> 

1500 IFC 1.E-5-ABSC l.-X2/X2NEW>> 11 

1510 I FC 1 • E- 5- ABSC 1-Y2/Y2NEW) ) 11 

1 520 X2«X2NEW+»2*CX20-Xa'JEW) J Y2sY2NEW+.52*CY20-Y2NEW) 
1 530 Z*Z*H 

1540 CALL0UTPCX2*Y2*Z> 

1 550 IFC ABSC ST0P-ZO5-ABSCZ+. 00 l*H-Z0>) BEGIN* BEGIN 
1560 X0»XllYO»YliXlaX2iY l laY2 
1570 G0T0 55 

1580 11 X2»X2NEWJ Y2*Y2NEWj GOTO 10 
1590 SUBROUTINE TITLE 



PMTEXP C0NTINUEO 



t600 PRINT* PRESSURE V0LUM E TEMPERATURE HEAT FL0W** 

U10 PRINT' C ATM* ) (LITERS) (DEC KELVIN) (CAL0RIES)" 

1680 RETURN 

J 630 SUBR0UTINELINE(P* V* T* Q) 
1640 C0MM0N FF(4) 

1650 PRINT* P*FF( 1), V*FF( 2) * T* Q*FF( 4) 
16$0 RETURN 
1670 SDATA 

1630 1 4 242484 1 

1690 9.86927E-6 .999972 1. 2.3912E-4 
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^Language: ""FORTRAN)* 

Comments 

5 is a vector containing the current- values of P,V f T and Q to which it 
is equated by the EQUIVALENCE statement. H is a set of wpight factors 
used in the integration routine and FF is a set of conversion factors 
to allow for the fddt that MKS units are used internally and 
atmospheres, liters, and calories are used for input and output. 

multiplying by FF(1) converts nt/m 2 to atmospheres 

multiplying by FF (2) converts m 3 /kgmole to li ters/gm mole 

multiplying by FF (3) converts from °K to °K 

multiplying by FF (H) converts from joules to calories. 

F is the gas constant. A 5 and F5 are the Van der Waals constants, a 
and b, and may be easily changed here. 

The present code for substance number is: 1 = ideal gas; 2 = van der 
waals gas; 3 = ideal Debye solid. If the substance number is not 
positive, the program stops. 

LFTTR convents letters to numbers such that: 

P becomes 1 
V becomes 2 
T becomes 3 
Q becomes 4 

The numbers are then used as identifying parameters. 

A soluticn of the equation of state for V(P,T). Other cas*s may be 
entered here if desired. The iterative solution for the Van der Waals 
r*z~: is suggested by Tribus. No attempt has been made to select an 
appropriate root when below the critical point. The convergence 
sometimes fails in this region., Therefore r to study the Van der Waals 
equation in this region, one should specify the volume. 

The volume of the solid is calculated from 

d (In V) = - kappa * dP + beta * dT 

kappa = (680 + .275 * T) 10- 9 /atm 

beta = (4 1. 5 + .023 *T) * 10-*/°C at P = 1 atm. 

(These values are from Allis £ Herlin p. 32, for copper.) The equation 
of state follows after converting to MKS values and integrating from 
the assumed initial state (1 atm., 7,1*10- 3 m 3 , 300°K) first along P=1 
atm. until T - Tf, then along that isotherm until P = P f . Clearly this 
stretches the data beyond its validity, but the result may be used for 
instructional purposes. 

The pressure solutions. 

The temperature solutions. 

This routine supplies Cp, Cv f beta and kappa. Beta (the coefficient of 
expansion) and kappa (the compressibility) are obtained from the 
eguations of state. For the ideal gas, C p and C v are tho'se for a 
rnonotonic gas and independent of temperature. For the Van der Waals 
gas, Cp is as given by combining the value for an ideal diatomic gas 
(Tribus, p. 177ff, h'2 "values) with the Van der Waals corrections given 
by Crawford, p. 338f£ . The hyperbolic term is completely negligible 
below 50°K. The Debye function expansions used for the solid are as 
given by Zharkhov and Kalinin, p. 2M . Numerical tests show this 
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calculation |to give nearly four-figure accuracy in 
is near the Debye temperature). 
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This routine sets the derivative to be used in the Simpson rule 
in teg rati on. The user has specified which variable is to be const ant 
and which is the controlled variable. If only p, V, and T are 
involy.ed, the other variable is determined from the equation of state 
and the appropriate derivative for 0 may be expressed (by the usual 
methods) in terms of P, V, T, C , C v , beta and .kappa. Which derivative 
is required is determined by the value N2, set previously. 



This routine determines P, V, or T from the equation of stat<=», 
variable to be determined is specified by the value of L. 
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,Tf Q is one of the specified param 
variables must be determined by jn 
derivatives. They are determined 
particular case specified by N2, w 
the fixed and stepped variables, 
variable is P, V, or 7, respective 
For M2 = 10, 11, or 12, the fixed 
and Q is the stepped variable. In 
be determined by integration. F i 
first required integration and G i 
(The order is P, V, then T) . 



eters, then two of the P-V-T 
tegration. F and G are the required 
from the Td 5 equations for the 
hich carries information about boMi 
For N2 = 7, 8, or 9, the stepped 
ly , and the process is adia bat ic. 
variable is-P, V or ?, respectively, 
each case, two P-V-T variables must 
s the appropriate derivative for the 
s the derivative for the second. 



Since the variables X, y, and Z may be one of several arrangements of 
the set P,V,T, and Q, it is necessary to. sort. ' ! ->m into proper order 
before printing. The selection of this order is determined by the 
value of N2. 

This is the first part cf the input routine. Th° equation is read, 
ignoring the = sign, the proper subscript is determined by LETTR, 
and the conversion factors are used. 

If either I or' J was 4 (i.e., Q) , another initial par ameter ' is needed. 
Otherwise, N2 is set according to the following tabl° and the remaining 
P-V-T value is obtained from the eguation of state. 
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Prints a title, intiali7.es Q, and prints the initial line. 

This line initializes the variable used in the Simpson rule integration 
and also stores the starting value for later reference. 

L becomes the subscript for 'the P-V-T value which will be determined 
from the eguation of state. 

A nine-point Simpson rule integration (see Pennington, p. 1 9 3 f f ) for 
Q. S(J) is + he stepped variable and S(L) is determined from th° 
eguation of state. The other variable, S(I) is fi xed . The stopping 
criterion in line 1260 allows for both negative and positive steps and 
for e small round-off error in the value of VAB. 



1290 - 1310 
1320 
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1340 



This is similar to 1110 in case another P-V-T value was required. 

If J is not 4, then 1=4 here. The value of N2 is set and X0 and. Y0 

are the (ordered) variables to be determined by integration. 

If J is 4; similar to 1320. 



In it ia li ze 
li ne. 



the stepped variable for integration and print the initial 
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1350 - 1580 The integration of two simultaneous differential equations foi the two 

variables (P, V or T) not determined by the input, simultaneous 
equations are required since the slope at intermediate points alonq 
the path may depend upon the value of the other variable at those 
intermediate points. The method used here is taken from McCracken G 
\ Dorn, p. 343ff. 

. 1590 - 1620 Prints title, 

1630 - 1660 Prints a line of output data after converting to appropriate units. 

1670 - 1690 The conversion factors in 1670 may be changed if another system of 

units for 1-0 is desired. 
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